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P\) ! Abstract 

In the London limit of the Ginzburg-Landau theory (Abelian Higgs model), vortex dipoles 
(small vortex loops) are treated as a grand canonical ensemble in the dilute gas approximation. The 
summation over these objects with the most general rotation- and translation invariant measure 
of integration over their shapes leads to effective sine-Gordon theories of the dual fields. The 
representations of the partition functions of both grand canonical ensembles are derived in the 
form of the integrals over the vortex dipoles and the small vortex loops, respectively. By virtue 
of these representations, the bilocal correlator of the vortex dipoles (loops) is calculated in the 
Ph! low-energy limit. 

It is further demonstrated that once the vortex dipoles (loops) are considered as such an 
ensemble rather than individual ones, the London limit of the Ginzburg-Landau theory (Abelian 
Higgs model) with external monopoles is equivalent up to the leading order in the inverse UV cutoff 
to the compact QED in the corresponding dimension with the charge of Cooper pairs changed 
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due to the Debye screening. 
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1 Introduction 



It is commonly argued that Abrikosov vortices ||T[ (Nielsen-Olesen strings @) in the dual Ginzburg- 
Landau theory (Abelian Higgs model) reveal the properties similar to those of strings in 3D- and 
4D QCD, respectively. This observation is based on the so-called 't Hooft-Mandelstam scenario 
of confinement , according to which confinement in QCD can be macroscopically thought of as 
the dual Meissner effect. Once being put forward, such a correspondence then enabled one to de- 
velop various phenomenological models of confinement, e.g. the so-called dual QCD approach [|J. 
However, an extremely strong support has been given to the 't Hooft-Mandelstam scenario by the 
so-called Abelian projection method || (see e.g. Refs. || [?|, [8|, |9], [10], |TTJ, [12], [L3|, Q for recent devel- 
opments and Ref. [15| for a review). The main outcome of this method is that under certain quite 
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pleasurable assumptions (like the so-called Abelian dominance hypothesis 
can indeed be with a good accuracy viewed as the London limit of the dual Abelian Higgs model 
with external electrically charged particles. An independent support of this conclusion comes out 
from the evaluation of field correlators in the dual Abelian Higgs model and comparison of them 
with those in QCD, introduced within the so-called Stochastic Vacuum Model []T7] and measured in 
the lattice experiments in Ref. fL8|. This calculation has been performed in Refs. fl9|, ED, 2j|, and 



as a result a very good agreement between the two sets of correlators has been established Q. All 
that encourages one to work further on investigating the confining properties of the Abelian Higgs 
model, which is a good and simple example of a model exhibiting the property of confinement. 
This is the main motivation for the present research. 

In a variety of an existing literature on the Abelian Higgs model from the point of view of 
the Abelian projection method (see e.g. Refs. P2] , |2"3"| , |2"T]|), the interaction of the Nielsen-Olesen 
strings has not been studied. In this respect, string representations of various objects (like partition 
function or field correlators) obtained there are insensitive to the properties of the ensemble of 
strings as a whole and depend actually on a single string only Q. Apart from that, there exist some 
publications where vortices in the 2D- and 3D Ginzburg-Landau theory have been treated as a 

26| . The results of these investigations are as follows. In 



grand canonical ensemble of dipoles 
the 2D case 



25 1, the effective field theory, emerging after the summation over the vortex dipoles, 



is the sine-Gordon one with the action of the type (] 
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Here, fh stands for the gauge boson mass, and C, is the so-called fugacity of dimension (mass) 2 , 
which is the statistical weight (Boltzmann factor) of a single vortex dipole. As far as the 3D 
case [^(J is concerned, there the resulting effective field theory is again the Ginzburg-Landau type 
theory of magnetic Higgs field (albeit with an additional mass term of the dual vector field h), 
whose partition function has the form 



1 A similar calculation of field correlators in the Abelian-projected 5?7(3)-gluodynamics has been performed in 
Ref. @. 

2 An obvious generalization of the string representation for the partition function of the theory with a global 
[/(l)-symmetry to an ensemble of noninteracting strings has been performed in Ref. J24|. 

3 Throughout the present paper, we work in the Euclidean space-time. 
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Here, if^j, = c^/i*, — d u h^ is the field strength tensor of the dual field h, r\ is the v.e.v. of the 
original Higgs field, which describes electric Cooper pairs, whose charge q is the double electron 
one. In order to arrive at Eq. (0), one should sum up over the grand canonical ensemble of the 
vortex dipoles, specifying the path-integral measure to be the one of the gas with a short-range 
repulsion as follows: 
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Although both the free part of the world-line action standing in the exponent on the R.H.S. of 
Eq. (^D and the 5-type interaction are quite natural, it looks desirable to perform the summation 
over the grand canonical ensemble of the vortex dipoles with the most general rotation- and 
translation invariant integration measure without specification of its concrete form. Clearly, such 
an approach would then enable one to perform also the summation over the grand canonical 
ensemble of the small vortex loops, built out of the Nielsen-Olesen strings in the Abelian Higgs 
model. Notice that the summation over the vortex loops seems to be difficult to perform by 
a direct generalization of Eq. (|3]), since the string world-sheet is a 2D object, rather than a 
ID vortex line. It is our aim in the present paper to proceed with such a summation over the 
small vortex loops in the Abelian Higgs model. As far as the summation over the vortex dipoles 
in the Ginzburg-Landau theory is concerned, it turns out to be analogous to the summation 
over the grand canonical ensemble of magnetic loops emerging in the Abelian-projected SU(2)- 
gluodynamics |14]] . In both cases of the Ginzburg-Landau theory and Abelian Higgs model, the 
effective field theories resulting from the summation occur to be of the sine-Gordon type (|l|). They 
allow for the representations directly in terms of the integrals over the vortex dipoles (loops). As 
we shall see, besides the Biot-Savart type interaction of the objects under study, the effective 
actions, corresponding to these representations, contain a multivalued potential resulting from 
the cosine interaction. In the low-energy limit, the real branch of this potential takes a simple 
quadratic form, and the integration over the vortex dipoles (loops) becomes simply Gaussian. 
This will enable us to calculate exactly the bilocal correlation functions of the vortex dipoles and 
small vortex loops in this limit. After that, we shall investigate the effects brought about by 
such a way of summation over the grand canonical ensemble of the vortex dipoles (loops) to the 
Ginzburg-Landau theory (Abelian Higgs model) with external monopoles. In particular, it will 
be demonstrated that to the leading order in the inverse UV momentum cutoff, this method of 
summation leads to an equivalence of the so-extended Ginzburg-Landau theory (Abelian Higgs 
model) to the compact QED in the corresponding dimension. 



The paper is organized as follows. In the next Section, we start our analysis with the Ginzburg- 
Landau theory and then apply the so-developed techniques in Section 3 to the Abelian Higgs 
model. The main results of the paper are summarized in Conclusions. Finally in three Appendices, 
some technical details of the performed calculations are outlined. 

2 Grand Canonical Ensemble of the Vortex Dipoles in the 
Ginzburg-Landau Theory 

Let us consider the partition function of the 3D Ginzburg-Landau theory in the London limit, 
i.e., the limit of infinitely heavy Higgs field, 
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where from now on (both in the case of the Ginzburg-Landau theory and Abelian Higgs model) 

K[ |23) 0, gg, 0. Here, F^ = d^A v - d u A„ is the field 
and 9 = 6 |smg ' + 6 reg ~ stands for the phase of the Higgs 



we adopt the notations of Refs 

strength tensor of the vector potential A^, 



field, consisting of a singular and a regular parts. First of them describes a certain configuration 
of Abrikosov vortices (e.g., a single vortex) according to the equation 
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(5) 
f dx„(r)5(x — x(r)) is the ^-function defined w.r.t. the vortex line parametrized 



by the vector x(r), < r < 1. The vortex lines are closed, which is expressed by the equation 
dpS^ = 0. Notice that Eq. (|5|) is nothing else but the local form of the Stokes theorem. As far as 
the usual regular part of the phase of the Higgs field is concerned, it describes fluctuations above 
this vortex background. 

One can perform the so-called path-integral duality transformation (see the above mentioned 
Refs.), which eventually casts the partition function (U) into the form of the integrals over the 
dual magnetic field h and vortex lines. Referring the reader for details to the above cited papers, 
we present here the final result of this transformation, which has the form 



Z, 



3D 
G-L 



T>~x.{t)T>\\D^) exp 



d?x 



v Elv + \T2 K + ^ 



27ri/i M <^ 



The scalar field tp can be further eliminated by performing the gauge transformation, h^ 
hfj, — —d^ip, after which the partition function takes the form 
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Notice that all the manipulations with Eq. @, which will be performed below, should obviously 
be done at the saddle-point for the field h, where due to the closeness of vortices, the second Proca 
equation, d^h^ = 0, holds. 

One can now sum up over the vortex dipoles in Eq. (||) in the sense of Eq. @, which leads to 
Eq. (E3). Our aim here is, however, to perform such a summation not specifying the measure to the 



particular form of the R.H.S. of Eq. (|3|), but rather keeping it as the most general rotation- and 
translation invariant one. To start with, we shall consider partition function (|6]) as a contribution 
of N vortex dipoles to their grand canonical ensemble. This can be done by replacing 5^ in Eq. (||) 
by the following expression describing the density of the gas of the vortex dipoles 
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Here, n a 's stand for winding numbers, and we have decomposed the vector x a (r) as follows, 

1 
x a (r) = y a + z a (r), where y a = / dTx. a (r) denotes the position of the a-th vortex dipole. In what 


follows, we shall restrict ourselves to the minimal values of winding numbers, n a = ±1, which 

is just the essence of the dipole approximation. That is because the energy of a single vortex is 

known to be a quadratic function of the flux [|IJ , due to which the existence of a dipole made out 

of two vortices of a unit flux is more energetically favorable than the existence of one vortex of the 

double flux. Besides that, we shall work in the approximation of a dilute gas of the vortex dipoles. 

' 2 which we 



According to this approximation, characteristic sizes of the vortex dipoles, / dr^ 


shall denote by a, are much smaller than characteristic distances |y| between them, which we 

shall denote by L. In particular, this means that the vortex dipoles are short living objects. 

Within these two approximations, by substituting Eq. (|7|) into Eq. @ one can proceed with 

the summation over the grand canonical ensemble of the vortex dipoles. This procedure essentially 



parallels a similar one of Ref. [J 3], and its details are outlined in Appendix A. As a result, we 
arrive at the following expression for the grand canonical partition function 
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Here, ( oc e~ So is the so-called fugacity, which has the dimension (mass) 3 with So standing for the 
action of a single vortex dipole, and A is the UV momentum cutoff. Thus, the summation over 
the grand canonical ensemble of the dipoles, built out of Abrikosov vortices, with the most general 
form of the measure of integration over their shapes yields in the dilute gas approximation the 
effective sine-Gordon theory (H). In particular, this way of treating the gas of the vortex dipoles 
leads to increasing the mass of the dual field h. Namely, expanding the cosine in Eq. (§) we get 
the square of the full mass, M 2 = m 2 + m 2 D = Q 2 1] 2 , where m = qrj is the usual mass of the 
dual field (equal to the mass of the gauge boson), and mu = fV^C is the additional contribution 

coming from the Debye screening. We have also introduced the full electric charge Q = Jq 2 + j^. 
Our next aim is to derive the representation of the partition function (^) directly in the form 
of an integral over the vortex dipoles. This can be done by making use of the following equality 
(valid at the saddle-point of the field h, corresponding to the action standing in the exponent on 
the R.H.S. of Eq. (|)), 
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Substituting it into Eq. (g), one can straightforwardly resolve the resulting saddle-point equation 
for the field h, -r^ sin (M) = —^j^jfi, which yields the desired representation in terms of the vortex 
dipoles 
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Here, the complex- valued potential of the vortex dipoles reads (c/. Ref. ||14 
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The obtained representation @ can now be applied to the calculation of correlators of the 
vortex dipoles. Indeed, it is possible to demonstrate that if we introduce into Eq. (|6|) (with 5^ 
a unity of the form 



replaced by 5^ as 
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and integrate out all the fields except j, the result will coincide with Eq. @. This is the reason why 
the correlators of j's are nothing else, but the correlators of the vortex dipoles. Such correlators 
can be calculated in the low-energy limit, i.e., when A|j| <C C an d one considers in Eq. fllQ|) the 
stable minimum of the real branch of the potential by extracting from the whole sum the term 
with n = 0. In this case, the bilocal correlator of the vortex dipoles reads 
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This example illustrates how the vortex dipoles in the grand canonical ensemble are correlated to 
each other. In particular, we see that their correlator decreases exactly according to the Yukawa 
law with the screening provided by the full mass M. 

Let us now see what will be the consequences of accounting for interaction of the vortex 
dipoles in the grand canonical ensemble if we introduce into the Ginzburg-Landau theory external 
monopoles. To start with, we shall introduce the monopoles into the Ginzburg-Landau theory 
with noninteracting Abrikosov vortices. This can be done by replacing the field strength tensor 
F^v in Eq. (|j) by F^ u + F^f,, where the monopole field strength tensor F^ obeys the equation 
\^^u\d^F^l = 2np with p standing for the density of monopoles. The path-integral duality 
transformation of the so-extended theory (Q) has been performed in Ref. |13|| and effectively 
results to adding to the Lagrangian standing in the exponent on the R.H.S. of Eq. (01) the term 
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We are now in the position to investigate what will be the effect of summation over the grand 
canonical ensemble of interacting vortex dipoles to the Ginzburg-Landau theory with external 



, i.e., keep in 



monopoles. If we restrict ourselves to the leading term in the -^--expansion of Eq. 

the expansion of the cosine only the term quadratic in h, than the integration over h is Gaussian 

and the result has the form 
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The details of a derivation of Eq. fll3|) are presented in Appendix B. This equation means that in 
the physical limit, when the UV momentum cutoff infinitely increases, the partition function of 
the Ginzburg-Landau theory with external monopoles, where the vortex dipoles are summed up 
in the sense of the grand canonical ensemble, is equivalent to the statistical weight of 3D compact 
QED with the charge q of Cooper pairs replaced by the full one Q. 

Finally, it is worth noting that accounting for all the terms in the -^-expansion of Eq. (§), 
rather than only for the leading one, leads simply to the following substitution in Eq. (W) 
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Grand Canonical Ensemble of the Small Vortex Loops 
in the Abelian Higgs Model 



In the present Section, we shall extend the analysis of the grand canonical ensemble of the vortex 
dipoles to the 4D case of the small vortex loops, built out of Nielsen- Olesen strings, in the London 
limit of the Abelian Higgs model. The partition function under study is given by Eq. (|4|) with 
the replacement d 3 x — > d 4 x. Notice that in what follows, we shall mark the 4D quantities (e.g., 
charge of Cooper pairs, v.e.v. of the Higgs field, etc.) by prime in order to distinguish them from 
the corresponding 3D ones. Equation (|5]) goes over into [^4], ^2], ||, |H], ^3|, [L3|] 
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J da pu (x(^))5(x — x(£)) is the so-called vorticity tensor current defined at the 

E 

string world-sheet S, parametrized by the vector x M (£). Here, £ = (£ 1 ,£ 2 ) G [0, 1] x [0, 1] denotes 
the 2D coordinate, and da fiu (x(^)) stands for the infinitesimal world-sheet element. Due to the 
closeness of strings, the vorticity tensor current is conserved, i.e., d^T,^ = 0. 

The path-integral duality transformation of the Abelian Higgs model accounting for noninter- 
acting Nielsen- Olesen strings has been performed in the above mentioned papers, and the result 
reads 
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Here, H^x = d p h uX + dxh^ u + d u h Xp , stands for the field strength tensor of the massive antisym- 
metric tensor field h^ v , usually referred to as the Kalb-Ramond field [27j. Analogously to the 



3D case, the vector field B^ can be eliminated by performing the hypergauge transformation, 
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In order to proceed from the individual string to the gas of the small vortex loops, one should 
analogously to the 3D case substitute for E MJy in Eq. (|T4]) the following expression 
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After that, the summation over the grand canonical ensemble of the vortex loops is straightfor- 
ward. Referring the reader for the details to Appendix C, we shall present here the result of this 
procedure, which has the form 
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Here, \h^ u \ = Jh 2 , A' is a new UV momentum cutoff, and the fugacity £' (Boltzmann factor 

of a single vortex loop) has now the dimension (mass) 4 . The square of the full mass of the field 
h^ following from Eq. ([[(]) reads M' 2 = m' 2 + m'^ = Q' 2 r]' 2 . Here, m' = q'rj' is the usual Higgs 

contribution, m' D = -^ is the Debye contribution, and Q' = J q' 2 + pj is the full electric 
charge. 

The representation of the partition function ( |T6"D in terms of the vortex loops can be obtained 
by virtue of the following equality valid at the saddle-point of Eq. (|T4"D (with S Mi/ 
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where Ki stands for the modified Bessel function. Substituting this equality into Eq. (P^|), we can 
integrate the field h^ v out, which yields the desired representation for the partition function (|i~6"D, 
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where the effective potential V is given by Eq. (|II]) with d 3 x — > d A x. 

Similarly to the 3D case, correlation functions of S^'s, calculated by virtue of the partition 
function (0), are nothing else, but the correlation functions of the small vortex loops in the gas. 
This can be seen by mentioning that if we insert into Eq. fll~4] ) with E Ml , replaced by E^f the 
following unity 
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and integrate out all the fields except S^, the result will coincide with Eq. ([IT]). Such correlation 
functions of the small vortex loops can be most easily calculated in the low-energy limit, A' 2 \Sp U \ <C 
(', by considering the stable minimum of the real branch of the potential (|10|) , where it takes a 
simple parabolic form. In particular, the bilocal correlation function reads (cf. Eq. (|TT|)) 
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Thus, we see that if Nielsen- Olesen strings are considered not as individual ones, but as a gas of 
the small vortex loops, the interaction between which has the Yukawa form, the bilocal correlator 
of the vortex loops in the low-energy limit has also the Yukawa behaviour (albeit screened not by 
the mass of the gauge boson, but by the full mass M'). 

Let us now investigate grand canonical ensemble of the vortex loops in the presence of external 
monopoles. To start with, we first clarify what will be the modifications of the Abelian Higgs 
model with Nielsen-Olesen strings treated as individual ones due to external monopoles. The 
starting partition function has the form 
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Here, the monopole field strength tensor F^ obeys the equation d lx F t 
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with jp standing for the (conserved) monopole current. The path-integral duality transformation 
of the partition function (|T8| ) has been performed in Ref. |21| and leads to Eq. flI3|) with the 
substitution S Mly -> £ Mi , = E^ + t^F^' . 

It is now straightforward to see what will be the consequences of the summation over the grand 
canonical ensemble of the small vortex loops to the Abelian Higgs model with external monopoles. 
Namely, this summation leads to the following expression for the partition function: 
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In particular, in the physical limit A' — > oo, keeping in Eq. fl2"0|) only the leading term in the 
■^■-expansion, we obtain 
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The derivation of Eq. (|2ID is similar to that of Eq. (|TB|). The only technical detail necessary for it 



is the integration over the Kalb-Ramond field, which can be found in Refs. p8], ^1[ . 

Thus we see that similarly to the 3D case, in this physical limit, the partition function of the 
grand canonical ensemble of the small vortex loops extended by external monopoles is equivalent 
to the statistical weight of the 4D compact QED with the electric charge q' replaced by the full 
one Q'. 

Finally, accounting for all the orders in the l/A'-expansion in Eq. (|20"D is equivalent to the 
following substitution in Eq. flTT|): 

1 ~ 
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This means that external monopoles do not affect the Yukawa-type interaction of the vortex loops, 
but enter only their effective potential. 

4 Conclusions 

In the present paper, we have investigated grand canonical ensembles of the vortex dipoles and 
small vortex loops in the London limit of the Ginzburg-Landau theory and Abelian Higgs model, 
respectively. In the approximation where these objects form a dilute gas, the summation over 
them with the most general rotation- and translation invariant measure of integration over their 
shapes has been performed. The resulting effective theories turned out to have the form of the 
sine-Gordon theories due to the additional cosine interaction term of the dual field, emerging from 
the above mentioned summation. In the physical limit of the large UV momentum cutoff, one can 
keep in the expansion of this cosine interaction in powers of the inverse cutoff only the first term. 
This yields a certain (positive) correction to the mass of the dual field due to the Debye screening. 

After that, we have casted the obtained sine-Gordon theories, corresponding to the Ginzburg- 
Landau theory and Abelian Higgs model, into the forms of the representations in terms of the 
vortex dipoles and small vortex loops, respectively. The resulting effective actions turned out to 
have similar forms and consist of the Biot-Savart Yukawa type interaction of the objects under 
study and a certain multivalued effective potential. In the low-energy limit, this potential takes a 
simple quadratic form, and the effective actions become Gaussian. This enables one to calculate 
correlation functions of the vortex dipoles (loops) in this limit, and as an example, we have 
calculated the bilocal correlators. Those turned out to have the Yukawa form, where the screening 
is governed by the full mass of the dual field, resulting both from the Higgs and Debye effects. 

Then, we have addressed the problem of what will be the consequences of our approach to 
treatment the vortex dipoles (loops) as a grand canonical ensemble to the theories under study 
extended by external monopoles. In this way, we have demonstrated that to the leading order 
in the expansion in powers of the inverse momentum cutoff, summation over the grand canonical 
ensemble results to the complete equivalence of these theories to the compact QED in the corre- 
sponding dimension with the charge of Cooper pairs replaced by the full one (which accounts also 
for the Debye screening). This result differs from that, which one gets by considering Abrikosov 
vortices (Nielsen-Olesen strings) as individual ones. In that case, the equivalence to compact QED 
holds only in the limit of vanishing gauge boson mass. As far as the effects brought about by 
external monopoles to all the orders of the expansion in the inverse powers of the UV momentum 
cutoff are concerned, those result into certain shifts of the arguments of the effective potentials, 
but do not affect the Biot-Savart interaction. 
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It now looks attractive to investigate the grand canonical ensembles of the vortex dipoles 
(loops) emerging in the effective Abelian-projected theories, corresponding to the original SU(N), 
N > 2, Yang-Mills theories, since the dual Abelian Higgs model can be under some assumptions 
considered as such an effective theory for the case N = 2. Recently, some progress in this direction 



has been achieved in Ref. 30 
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Appendix A. Summation over the Grand Canonical Ensem- 
ble of the Vortex Dipoles 

In this Appendix, we shall present some details of a derivation of Eq. (^ of the main text. Sub- 
stituting Eq. (0) into Eq. (j^j, one can perform the summation over the grand canonical ensemble 
of the vortex dipoles as follows: 

oo /-N / N . „ \ ( N 



1 + E W\ ( ft / rf V / Vz*n M ) £ exp J 2m ±nj dz^h, (x a ) 

> ■ iv - \i=l J J J n a =±l { a=l 

1 + E ^§r {/ d3 y J V M*\ cos (277 j dz,h, 



-/i"7i(, X / 



JV=1 

Here, we have introduced the so-called fugacity (, which is proportional to the statistical weight of 
a single vortex dipole and has the dimension (mass) 3 . Next, /x[z] stands for a certain rotation- and 
translation invariant measure of integration over the shapes of dipoles. Let us now employ for the 
grand canonical ensemble of the vortex dipoles under consideration the dilute gas approximation 
described after Eq. (|7]). One can then expand /i M (x) up to the first order in a/L as follows: 



M x ) = My) + L- l z u n u h^(y) + O ({^)\ (Al) 



where n u = y u /\y\, and we have substituted n u /L for the derivative d/dy u . By making use of this 
expansion, we obtain 

Vz/j,[z] cos ( 2ix j> dz M M x )J - / "^z/ifz] cos ( -r n uh fl (y)P lll y[z 

(-1)" /2tt\ 
„_o (2")« V L J 



00 (—1)™ /2ir\ 2n f 

E T^r (~l) n ^ h ^ y } ' ■■ n ^Mn(y)y v M z } p ^A z ] ■ ■■p^^M 



~ (-1)" /2vra 2 V n 

E T^yr I — ^~ J n "i V ( y ) ' ' ' n ^ v (y ) x 
i 

t/j 1 I/ 1 ,/i 2 l'2 ' ' ' *-fJ.2n-lV2n-l, A»2n^2n "T pemiUtatlOnS 



(2n-l)!! 



ff 



Here, -P^z] = § dz p z u is the so-called tensor area, and l^xp = \ (SpxSvp — 8p P 8u\)- Notice that 
the sum in square brackets on the R.H.S. of the last equality contains (2n — 1)!! terms, which is the 
reason for extracting explicitly this normalization factor. The general form of the tensor structure 
standing in these brackets is due to the rotation- and translation invariance of the measure /x[z]. 
In the further contraction of indices, it is worth noting that within the dilute gas approximation, 
n fihfj,(y) = 0. Indeed, as we have already seen, this approximation allows for the substitution 
n M — > L-t^-. On the other hand, the divergency of the field h vanishes according to the second 
Proca equation, valid at the saddle point of the field h corresponding to Eq. (|6|). This completes 
the proof of our statement. 

Finally, introducing the UV momentum cutoff as A = -M 2 (^> a^ 1 ), we obtain 



/ Dz/i[z] cos ( 2tt j> dz t 



A,fx)^cos^ |h(y ^ 



wv ~ v — ^ A J 

which leads to Eq. @. 

Appendix B. Details of a Derivation of Eq. (J3) 

In the present Appendix, some details of a derivation of Eq. fll3|) will be presented. Let us start 
with the following expression: 

J S J Ph exp { - / A [J^ + f h> + '-H , A / ^ 

which (up to an inessential constant factor, which can be referred to the integration measure) 
is the leading term in the ^--expansion of the partition function (|8]) including external monopole 
part (|i~2"P . Gaussian integration over the field h yields 



J = exp < — — / d 3 xd 3 y 



e -M|x-y| 

47r|x — y| 



d 2 f ,o p(z) r ,o p(u) Aw 2 . . , . 



AdxndyaJ |x — z|j |y — u| M 2 



(B.l) 
where the normalization factor was assumed to be included into the measure T>h. 

Let us consider the first term in square brackets on the R.H.S. of Eq. (B.l) together with the 
factor in front of these brackets. By performing the partial integration, one can cast the derivative 
d/di/fj, to the factor in front of the brackets, after which due to the translation invariance of this 
factor it can be replaced by —d/dx p . This derivative can be casted back by doing one more partial 
integration, which yields for this group of terms 



- / d 3 xd 3 y- 
AJ y < 



Air\y:-y\dx^dy^J |x-z|7 |y - u| 



1 fd 3 xd 3 y^ r A x fd 3 z-^- J d 3 u-^- = 7r l ' d 3 xd 3 yd 3 up{*) 

A J Air x — v J x — z \J v — u J 



y| J |x — z\ J |y — u| J 47r|x — y| |y — u| 

Substituting this expression into Eq. (B.l), we obtain 
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J = exp 
The integral 



TIT] 



d 3 xd 3 y 



d 3 y- t 



e -M|x-y| 

47r|x — y| 



-M|x-y| 



,o p(x)p(u) 47T 

^ I | + ^# 

|y — u M" 1 



^7 



-A/|y| 



(B.2) 



;b.3) 



|x-y||y-u| J |y||w-y| 

where w = x — u, can be calculated by the two alternative ways. One of them is to divide 



the integration region over |y| into two parts, [0, |w|], [|w|,+oo) and expand 



i 



|w-y] 



in Legendre 



polynomials P n 's on both of them. Then, the integration over the azimuthal angle singles out 

1 
from the whole series only the zeroth term, / P n (cos8)dcos8 = 25 n Q. After that, the integration 

-l 

over |y| at both intervals is straightforward and yields for the integral (B.3) the following result: 

47T 



M 2 lw| 



1 



-M|w| 



[BA) 



Another way to calculate the integral (B.3) is to use the relations 



d 3 p 



^py 



1 e 



-A/|y| 



(2tt) 3 p 2 + M 2 4vr |y| 
by virtue of which we have 



and 



d 3 q e 



iq(w-y) 



(2tt) 3 q 2 



d 3 pd 3 q e ipy+iq(w ~ y) 2 

" 167T 



(2tt) 6 (p 2 + M 2 )q 2 

-oo +oo 



d 3 p 



IQti 



I = 16tt 2 f d 3 y f 

d3p J da J ^ e ^pw- Q p 2 -/3(p 2 +M 2 ) =2 ^ J da J dp 



47r|w — y| 



«pw 



(2tt) 3 p 2 (p 2 + M 2 ) 

+oo +oo -)3M 2 - 



4(a+/3) 



It is now suitable to introduce new integration variables as a = bt, (3 = 6(1 — t); b 6 [0, +oo) 
t G [0, 1], which yields 



1 = 2^ f ^ dte-^~ bM2 ^ = 2 -^ I 
J \/h J M l J 



-j-oo 

2Jri f db _w£ , 
— — / — e 4t 1 

M 2 J b l 

o 



-bM' 



' 

The remaining integral is already straightforward to evaluate by changing the variable & — » | and 
making use of the formula 



db 
Vb 



e Ab -T 



T 



-2VAB 



valid for positive A and B. The result coincides with Eq. (B.4). 

Finally, substituting Eq. (B.4) into Eq. (B.2) and recalling that M = Qt], we arrive at Eq. (|T3J) 
of the main text. 
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Appendix C. Summation over the Grand Canonical Ensem- 
ble of the Small Vortex Loops 

In the present Appendix, we shall outline some steps of a derivation of Eq. (]T^|). Let us first 
consider the infinitesimal world-sheet element of the a-th vortex loop, which has the form 

*V (x a (0) = e a %d a x;(0)(dpx:(0)d 2 ^ 

where a, (3 — 1,2. Analogously to the 3D case, it is reasonable to introduce the center-of-mass 
coordinate (position) of the world-sheet y° = J ' d?£Xp(£). The full world-sheet coordinate can be 
respectively decomposed as follows x°(£) = y^ + <z"(£) with the vector z?(£) describing the shape 
of the a-th vortex loop world-sheet. Then, substituting into Eq. (|14D instead of E^ the vorticity 
tensor current of the vortex loop gas fll"5|), we can perform the summation over the grand canonical 
ensemble of the vortex loops as follows: 

H-ECfn/^V/^^VMl E exp{m£n a /aV(^(0)M* a (0)) = 

JV=1 iV ' \i=l J J / n =±l { a=l J ) 

= 1 + E ^^ {J d 4 y J Vz p (OAz] cos (^t J da,MO)h,u(x(0))} N - 

Here, /jf is a certain rotation- and translation invariant measure of integration over shapes of the 
world-sheets of the vortex loops . Employing now the dilute gas approximation, we can expand 
h,a, up to the first order in a'/L' (cf. Eq. (A.l)), where a' stands for the typical value of |z a (£)|'s 
(vortex loops sizes), which are much smaller than the typical value V of ||/ a |'s (distances between 
vortex loops) []. This yields 

Vz p (£)ii'[z] cos U j do>(z(£))V( x (0)J - j Vz pH'\A cos (jjn x h^ v (y)Pl vX [z 



00 ( — 1)™ / IT \ 2n f 

= E T^yr [jj) n Xl K lVl {y) ■ ■■n X2n h ll2nV2n (y) j Vz p ^\z\P'^ VxM \z\ ■ ■ ■ P^^xjz], (CI) 

where n x = y\/\y\ and P'^^z] = J dcr pv (z(^))z x (^). Due to the rotation- and translation invari- 
ance of the measure \J [z] , the last average has the form 

J Vz^'[z]Pl iUiM [z\ ■ ■■P' ll2nU2nM M = 



(a' 3 ) 2n 



L Ml t/ l.A t 2l / 2 AiA2 ' ' ' i A»2n-1^2n-l,A t 2n^2. 



^A 2n _ 1 A 2n + permutations 



(2n- 1) 
Substituting this expression into Eq. (C.l), we finally obtain 

4 Similarly to the 3D case, the approximation a' <C L' means that the vortex loops are short living objects. 
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oo / i\n /__/3^ 2n 



ir/WY II. /- M*> nV(»)l 



Jvz p (Ov'[z}cos (jr J da, u (z(0)KMO)) * £ 7^1" h^T IMv)!* 1 = «» 



n f (2n)! V # / ' " n V A' 2 



where we have introduced a new UV momentum cutoff A' = \/^m, which is much larger than 
1/a'. This yields the desired Eq. (|16|). 



wa" 
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